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In this article,, a semi -analytical theory of general planetary 
perturbations is developed, which is somewhat akin to Hill's 
theory. In both methods the first order perturbations coincide, 
but the theories of perturbations of higher orders are different. 

The inconveniences of Hill's method, namely: the triple integral 

in the perturbations of the radius vector and the redundant constant 
of integration, do not appear here. The short and the long period 
terms containing the squares of the small divisors are localized 
and combined together. The existence of such a direct way of 
separating these important terms from the remaining perturbations 
constitutes a significant characteristic of a planetary theory. 

The form of decomposition of perturbations as used in this article 
leads to a system of differential equations easily integrable by 
Hill's procedure and to a symmetrical scheme for the computation of 
perturbations of higher orders . 
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Notations 


m - the mass of the disturbed planet. The mass of the sun is taken as unity, 
m' - the mass of the disturbing planet, 
k - the Gaussian constant. 
fi 2 - k 2 ( 1 + m) . 

7 - the undisturbed position-vector of the planet m. 

° I r | . 

- the unit vector normal to the undisturbed orbit plane of the planet m. 

J - the angular velocity of rotation of the frame (r, Rx 7, r). 
r° - the unit vector in the direction of r. 
r' - the undisturbed position vector of the planet m\ 

Sr - the perturbation in the position vector of the planet m. 
r + Sr - the disturbed position vector of the planet m. 

Sr' - the perturbations in the position vector of the planet m'. 

7' + S7' - the disturbed position vector of the planet m\ 
v - the undisturbed velocity of the planet m. 

I - the undisturbed mean anomaly of the planet m. 
e - the undisturbed eccentricity of the planet m . 
n = — — the undisturbed mean motion of the planet m. 

a 3/2 

a - the undisturbed semi- major axis of the planet m. 
p = a (l - e 2 ) . 

f - the undisturbed true anomaly of the planet m . 

V - the del-operator with respect to r. 

V' - the del-operator with respect to r'. 

E - the base of natural logarithms. 

~p - 7 ' - 7 . 

7 k - the perturbations of the k th order in the position vector of the planet m. 
7 k * - the perturbations of the k th order in the position vector of the planet m\ 

A = • 

o(r , 7 # ) = - the main part of the disturbing function. 

I - the idemfactor. 
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Introduction 

In this article a numerical theory of general planetary perturbations is developed. The per- 
turbations are obtained in the standard form of series containing the periodic, the secular and 
the mixed terms. The coefficients of terms are numerical; they are obtained by double harmonic 
analysis as applied to the force components, and subsequent integration. The numerical theory 
of perturbations in the coordinates escapes the inconveniences of the numerical theory in the 
elements: for example, since the eccentricity and the sine of the inclination do not appear as 
divisors in the differential equations, no numerical difficulties arise in the case of nearly circular 
orbits. 

The method presented here is somewhat aton to HilTs method (1874), at least where the first 
order perturbations are concerned: in both methods the first order perturbations coincide. The 

theories of perturbations of higher orders in both methods are different. 

We determine the perturbations in rectangular coordinates directly without using the per- 
turbations in polar coordinates as an intermediary means, as is done in Hill's method. Some other 
inconveniences are removed in the theory presented here: the triple integral and the seventh 
constant of integration so peculiar to Hill's method do not appear; the components of perturba- 
tions along the radius-vector are determined in a more direct manner; and the difficulties as- 
sociated with determining the redundant constant of integration vanish. 

In Hill's method the undisturbed true anomaly is taken as the independent variable. Such a 
choice causes numerous inconveniences if the perturbations in the motion of the disturbing body 
are also to be taken into account. For this reason the universal variable, time, is used in our 
exposition. 

Since the advent of electronic machines, general perturbations theories can successfully 
compete with numerical integration procedures. The theory of Mars in Hansen's coordinates 
recently developed by Clemence (1949, 1961) brilliantly confirms this statement. The reference 
ellipse representing the undisturbed motion can be chosen in a variety of ways. The only restric- 
tion imposed is that the difference between the disturbed and undisturbed motions be small, of 
the order of the perturbations. 
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In the article by Musen and Carpenter (1953) a decomposition of the perturbations in the posi- 
tion vector along 7, v, and R was suggested. In Lhe present work we suggest a decomposition 
along 7°, Rx 7°, and r. This form of decomposition leads to a more compact and more symmetri- 
cal scheme for developing perturbations of higher order than the author's previous scheme. 


Basic differential equations. 


The equation of the motion of the planet m as referred to the rotating undisturbed frame 7°, 
R x 7° , R has the form 


d 2 87 -» d87 — I - \ m 

~^T + ^ x "3F + ^ x (i/r x 8r) + 3t‘ x 8r 


Substituting 



^iVf)(7 + 87, 7' + 87') 


(1) 


P 


r 2 


dip' : 2uVp dr - 


into (1), introducing the differential operator 

D = 87 • V + 87' • V' , 


and taking the equation 


1 87 377 • 

^>7 = -7T + — 7 


Sr 


into account, we can re-write (1) in the form 


d 2 87 2fj. Vp - 

— + ~7T-* 


dSr 
“cfT + 


M 2 P 


R x 


(Rx Sr) 


where 


2m Vp 


dr - 

■at R * 




m 2 f 


f = v(e d - i - d) i 


+ VE° n(r, 7') . 


( 2 ) 

(3) 
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In applying the operators V or V' we consider every function as a function of r and ? ' only and the 
perturbations Sr and S?' are considered as relative constants. The operator e d accomplishes the 
development of (3) into a power series with respect to the perturbations. 

Decomposing Sr a\;ng the axes T°, Rx 7 °, r and putting 


ST - f r° + 7? E x r° + £ R 


we have in the relative motion 


“dF 


d£ din - 

= Ht r° + St R * 


r ° + cft'R • 


d 2 Sr _ d 2 if; d 2 rj 


dt 5 


dt 2 


r° + — r R x r° + 
dt 2 


dH 

dt 2 


Substituting (4) through (6) into (2), we obtain 


d 2 <? 2 /u.yp d-o ( p 2 \ , 2 uyp dr _ , _ 

771 -at ^ + 7ij^ + — dt 7 ? - 



d 2 77 _ 2/j.yp d£ 2fiyp dr 

7777 -rfr-vt - 


f - 



m 2 h , 



m 2 z , 


(•i) 


(5) 

( 6 ) 


(7) 

( 8 ) 
(9) 


where we put 


• F , 


H = R x r° • F , 

Z = R • F . 

We shall make use of the standard method of developing the perturbations into power series 
with respect to the disturbing mass and put 


87 = r i + r 2 + r 3 + 


7 k = 4 r ° + t >k r * r ° + k R 
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(10) 


£ = + ^2 + ^3 + ••• ’ 

1) = 7 ?! + V 2 + V 3 + ••• ' 

C = 5i + C 2 + C, + ••• • 

F = F, t F 2 + F 3 + • • • , 

where r K> f K , tj k , i K , ? K are of the order k in m' . 

In the previous article by Musen and Carpenter (1963), on the basis of the formula 

VE D = V + (r, -VV+r,' • V' V) 

+ (r 2 • VV + ? 2 ' • V' V) + 3 (r 1 • V + 7/ • V') 2 vj 
+ (r 3 • W+r 3 ' • V'V) + (rj • V+7,' • V')(r 2 • V+7 2 ' • V ) V 


+ F ( 7 i • V + ? i' • V ') 3V 


( 11 ) 


it was found that 


„ m [p r ' 

F * = Vn = — 


( 12 ) 


< 3 ' 7 i 15 7 ( ?,7 i) 2 


2 


2 r s 

+ r J • vvn + r j' • V' vn , 


3 , 


F 3 ° ~s ( ? - 7 i 7 2 +?7 i • 7 2 + 7 , 1 ? 2 7 i) 


(13) 


15 __ _ 3 r i r i ’ r i 

~ 7 > rr i ■ rr 2 • r +_ 2 — 71 


15 77 ' ? i ? i • ? i 1S 7 i( ? - ? i) 2 


2 


35 'P-7,)* 
2 


2 


(r 2 -vvn +7; • v' vn) 


+ j (Tj • v + 7/ • v') 2 vn 


(14) 
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By substituting 


?i = r° + ^Rx r° + CjR 


and 


vvn = v 


1 + m ^ p 3 


T ' \ _ m' / _I_ 3pp \ 

‘ r' 3 J " 1 + E \ + ' 


P 5 / 


v' vn = V' 


1 + m ^ p 3 r '3 


m' f I 


1 + m \p* 


3pp 

nS 


1 + m V r ' 3 


into (13), we deduce a compact expression for f 2 : 


F 2 = [(' 3 ^i 2 + 4 7 7i 2 + 4 ^ 1 2 ) r° +37 ) l C 1 Rx r° + 3i 1 s 1 Rj 


Pi 3 pp • Pj ^ {tJ 3r ' r ' * r/ 

'3 r * 5 


1 **’ m i \ p$ j \ r 


Substituting (10) and (11) into (7)- (9) we obtain the equations 


d2 _ 2jj. -fp d ^K _ 2 /_p_ 2 \ dr 

at HE ^ y r 4 r 3 ) r 3 cTF 


, 2 S 


d 2 


2/x Vp 2/J.yp dr 

+ 1 -nr* : — tit 


dt 2 ' r 2 "3t r 3 <37 ~ ^ 


.r« r 3 / 




M 2 H k , 


d 2 C *■ 

b K 


dt 2 + 7 T ^ 


where 


* = 1, 2, 3, ••• . 


H* = R x r° ’ F K , 

Zk = R • F k ■ 


which are of the same form as (7)- (9). Thus, our problem now is to integrate the variational equa- 
tions of the form (7)- (9) with the right sides known. 
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Integration procedure. 


We shall make use of the substitution 


dr j 

£ - ur + dt |( w ” 2u ) 

(15) 

u yp f 

V = “f— J (w - 2u) dt 

(16) 

t = c . 

(17) 

which reduces (7)- (9) to the form integrable by Hill's procedure. Substituting 
(7), and making use of the relation 

(15) and (16) into 

ill . 2 /p i\ 
dt* ■ * [T> ~ 7*) 


We obtain 


d 2 u fi 2 1 dr dw fi 2 S 

+ 7 ? ( u - 2 w) + T-at dt = ~r • 

(18) 

From (15) and (16), we have 


1 / dr/ dr\ 2£ 

W ' mVp V dt ~ v dt j + r • 

(19) 

Differentiating (19) and taking (8) into consideration we obtain 


dw fir H 

dt " Vp 

(20) 

and 


* ■ K - * I VP dt ' 

(21) 


where the integral sign represents the integral obtained in a formal manner; K 3 is the constant 
of integration. 

Taking (20), (21) and the equation 


dr 

dt 


fie sin f 

Vp 
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into account, we obtain from (18): 


H 2 M 2 n 2 < r p «• i n f \ O C , . 3 r •* ; 

- 77 ^ fU - 2 a J r — fu - -Uv j - r 5 Q A / + — . — Gt . 

at" v •*' r J ^ r \ P y r - j yp 

The last equation can be integrated by using Hill's procedure and we have 


r r 

u - 2 X 3 + K ; 7 cos f + X 2 7 sir. f 


. \Ul- 




e s i n f 


r P 


H rr s 


in(f-f) 


dt 


sin (f-f)dt 


dt , 


(22) 


where k 1 and x 2 are constants of integration and r, f are considered as temporary constants; after 
the integration is completed they are replaced by rand f . 

The double integral in (22) can be simplified through integration by parts, and we obtain 


| a sin (f - f) dt | 


VrH 


dt 


2m rH 


yp 


dt 


2/i rrH 
P/P 


cos + e cos f J dt 


Taking this last relation into consideration we deduce from (22) after some easy transformations 

u = 2 K 3 + Kj - cos f + K 2 7 sin f + A , 

in which 


A = IMF • r t NF • Rx r dt 


where 


an r /— \ 

M = . ' ' rXsin(f-f) 




N 


(l 1 e 2 ) 


3/2 a 


1 T 1 

+ 2 e cos 1 " 5 e cos 


2 cos 


(f-2f) 
(f- f ) + 2 


The expressions M and n remain the same for the perturbations of all orders. They are to be de- 
veloped into a double Fourier series with respect to the mean anomaly l and with respect to the 
auxiliary mean anomaly l , associated with the auxiliary true anomaly f . 
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After the integration is performed l is replaced by l . We have 


“ cos fd£ 


3 

2 ent 


1 / r r 

2 K 1 ' e a sin f 




,2 a* 


sin f 


r 

T sin idl 

£i 



\ 

c cos 2 i } 

) 


(25) 

(20) 


Putting 


s n l (w - 2u) dt 


and taking (21), (24), (25) and (26) into account, we deduce 


s = 


3 K J. £l , , r 

n nt + n (3ent - -y/l - e a sin f - 


1 r- 


* 


,2 a 


sin f 


K 2 r 2 


“ e 2 — — (2 cos f + e cos 2 f ) + K 4 + B 


and, after some easy transformations, 


S " n (" K 3 + eK i) nt * 


r 2 


n yi - e 2 3 


2 [2 sin f + ~2 e sin 2f 


K2> / i T ^ K 2 r 2 f 1 1 

n In' “I 2 cos f + 2 e + j e cos 2f ) + K 4 + B 


(27) 


where 


narH 




-dt 


_ ^ 2 


1 “ e 


2A dt 


The formulas for the computation of 6 , tj and £ become 


(28) 


^ nae sin f 

f - ru + ' S 




V - 


YE 


■s , 


4 ” K s ~ cos f + K 6 “ sin f + C 


(29) 

(30) 

(31) 
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where C is the standard expression 


c 


r rr 



dt 




exrressior 


oi 


* A 


;'.iv _ i- 


are 


w- vvv^y 


C C * U «.C Ct 


math; 


is made. 


dt 


r 2 df 

® 2 n yi - e 2 


Determination of constants of integration . 

We consider here the determination of constants of integration for the case when the elements 
are osculating at the epoch t = 0 and for the case when they are mean. If the elements are os- 
culating, then we have 


(Sr), = 0 , (If), = 0 . (33) 

where the zero- subscript designates the value of the expression at the epoch. From (33) we 
deduce 


u o = ° . 


W 0 = 0 , 


l—) 

Wo 


= 0 


Taking (21) into account, we obtain 


$k 3 


naH 




-dt 


Differentiating (24) and taking the equations 


dr 

dl a cos f - 


sin f 


f 


d _r 
dl a 


sin f = + 


cos r -r e 




s 0 = 0 . (34) 


(35) 



into account, we deduce from (34): 


+ K. — cos f fl + K 2 — sin f 0 


••o —3 


sin * , 


cos f Q + e 


where wo put 


From these last equations we obtain 


cos f Q + e 


- q 2 ( A o + 2K 3 ) + " -- sin 


X * e 




sin f n 


(\ + 2k 3 




^0 r o 

T cos f 0 


In a similar way we deduce 


cos f Q + e 


C C ° r ° • r 

1 - e 2 C o T ~ = T Sln f o 




sin f 


o C o r o 

_ C o * r. ? T cos f o 


1 - e 




and putting t = 0 in (27) we deduce the following value for K. : 


k 4 = + 


K i r o 2 


n yi- 


.2 a 


2 sin f Q + 2 e 0 sin 2f, 


K 2 


n a 2 (2 cos f 0 + 2 e + j e cos 2f 0 j - B 0 


The mean elements can be defined in several ways. We accept here the following definition: 
the elements are mean if 

(1) the perturbations of the true longitude \ with respect to the orbit-plane defined by these 
elements do not contain the terms of the form 


K 0 , Kt, K< c) cosZ and K< s > sin l; and 


(36) 
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( 2 ) the expression for the M third coordinate" £ does not contain the terms of the form 


K (c) cos l , sin l 


We have 



cos mf 


2 C 0 n * m + Cj"’" 1 cos l + C 2 n ’ m cos 21 + 



sin mf 


S J n,m sin l + S 2 n ’ m sin 2 Z + • • • 


where 


C. n - m = X.' 


XJJ* 1 


S n = X n,m - X *!•' 


and x. n ' m are Hansen's coefficients. 

Let us start with the deter m i n a tio n of constants of integration in the perturbations of the 
first order. Perturbations of the first order in the true longitude are given by the expression v/r, 
where 77 = 75 ,, and consequently the terms (36) must be absent in the expression 


nr yT"-e 


If?) (" K 3 +eK l) nt 


K i ( 1 . A 

— — - 2 sin f + 77 e sm 2f 

n yr^ \ 2 1 

, 1 1 \ 

+ pj ( 2 cos f + 2 e + 2 e cos 2f ) 




(37) 


We have to separate the terms of the form (36) in the development of a a /r 2 B: 



a o + /3 0 nt + ctj cos Z + p x sin l + • • • 


(37’) 
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The condition for absence of the constant term in (37') leads to the equation 


K 2 yT- 


In a similar way we deduce 


( + C 0 «-> e ♦ \ eC».») + !<V-«K 4 ♦ a 0 = 0 


, ...» tC| oj , 

+ 'Si C o' J ’° (-K 3 +eK .) + /*0 = 0 


= 0 , 


= 0 , 


Separating in C the terms with the argument l , 


C = Cj cos l + Sj sin l + 


we obtain, taking (31) into account. 


$k,cM ♦ c, = o , 


( 38 ) 


(39) 

(40) 

(41) 


(42) 


K # S, 11 + S, = 0 . (43) 

In the planetary case the solutions of equations (38)- (43) can be found without any difficulty, 
because the coefficients only of one unknown in each equation are not small. The coefficients C"'" 
can be computed either by using the classical analytical expressions or by Cayley's tables (1861), 
or by means of harmonic analysis. The latter procedure is preferable if the eccentricity is not 
very small. 

Determination of constants of integration in higher order perturbations requires some addi- 
tional considerations. Let r,, r 2 , • • • be the perturbations in the radius-vector r and \ lf X 2 , • •• 
be the perturbations in the true longitude k of the first, second, etc., orders. From 


^ i+ £ 2 + ...) 7 o + 


K + 7 2 + 


N _ (dr dr\ 

•) Rx r° = ^r, + K^J 


+ 


dr d 7 \ 1 / 2 d 2 r d 2 r 

r 2 d? + S J + 2 ^ r i 777 + 2r i JT5K + 




a 2 r 
<?\ 2 


)1 
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and substituting 


<?r 

JT 


we obtain 


consequently, 


dr — _ 

dK = R * r ■ 


d 2 7 


dr 2 


d 2 7 

TFSX 


Rx r° , 



= r ! ' = rX I 


= r a ' 2 k i r ’ ~ tX 2 + r i • 


or, taking (30) into account, 


V 2 - if, V 1 


r 


X 2 = (8,-tf.SO , 

where S } corresponds to the first and S 2 to the second order perturbations. As we see, before 
the determination of the constants of integration of the second order in the case of mean elements 
a correction term -fj Sj must be added to b 2 . For the perturbation of the third order a similar 
correction term will depend upon £ lt s, , £ 2 , s 2 . 

Conclusion 

A revival of the general interest in planetary theories can be observed in our time. Several 
scientific institutions are dedicating their time and efforts to the astronomical solution of the 
planetary problem. The results by Brouwer (1944), Gontkovskaya (1958) and Danby (1962) must 
especially be mentioned. A considerable amount of work on the theoretical exposition as well as 
on programming has also been done at the Theoretical Division of Goddard Space Flight Center. 

In the present article we suggest a new scheme which is convenient for computing the perturba- 
tions of the first as well as of higher orders. 

The determination of constants of integration in the case of both the osculating and the mean ele- 
ments is a straightforward process in the proposed scheme. An important feature of the scheme is 
that the squares of small divisors, as caused by the commensurability of mean motions, are intro- 
duced by integration of only one expression, namely w - 2u. The short and the long period terms 
containing the squares of small divisors constitute a significant part in the perturbations £,v, i. 
The existence of such a direct way of separating these terms from the remaining perturbations 


constitutes a significant part of a planetary theory. The development here is kept in the form 
which facilitates comparison with the results obtained on the basis of the classical form of Hill's 
theory, if necessary. 
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